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What can Gaia (with TMT) say about Sculptor’s Core? 
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Walker et al.’s Magellan/MMFS Survey survey identified 1355 red giant candidates in the dwarf 
spheroidal galaxy Sculptor. We find that the Gaia satellite will be able to measure the proper 
motions of 139 of these with a precision of between 13 and 20 km/s. Using a Jeans analysis 
and 5-parameter density model we show that this allows a determination of the mass within the 
deprojected half-light radius to within 16% and a measurement of the dark matter density exponent 
7 to within 0.68 within that radius. If, even at first light, the TMT observes Sculptor then the 
combined observations will improve the precision on these proper motions to about 5 km/s, about 
5 years earlier than would be possible without Gaia, further improving the precision of 7 to 0.27. 
Using a bimodal stellar population model for Sculptor the precision of 7 improves by about 30%. 
This suggests that Gaia (with TMT) is capable of excluding a cusped profile of the kind predicted 
by CDM simulations with 2a (4<r) of confidence. 


According to the standard cosmological model, 
most of the matter in the universe is cold dark mat¬ 
ter (CDM), which at long distances only interacts 
gravitationally. Direct and indirect dark matter de¬ 
tection experiments are a major industry and so far 
have rapidly excluded large regions of the CDM pa¬ 
rameter space, as has the Large Hadron Collider. 
However such tests can only confirm CDM, any such 
attempt to falsify this paradigm can be evaded by 
fine-tuning. On the other hand, CDM makes fal- 
sifiable predictions for the density profiles of suffi¬ 
ciently dark matter dominated systems. They must 
fall as 1 /7' 3 at large radii and as at least 1/r at small 
radii, a scaling known as a cusp. While there is lit¬ 
tle consensus over just how much dark matter dom¬ 
ination is sufficient, studies such as Ref. mm sug¬ 
gest that galaxies with stellar masses beneath about 
10 6 ’ 5 ^ 7 Mq should be cusped while Ref. [3 argues 
that energy conservation implies that baryonic ef¬ 
fects cannot remove cusps in the Milky Way’s dwarf 
spheroidal (dSph) satellites except under quite spe¬ 
cial circumstances. Thus the discovery that a dSph 
is cored, not cusped, would pose a strong challenge 
to CDM. 

In this letter, following a strategy similar to that 
of Ref. 0, we determine just when it will be known 
whether the dSph Sculptor, with a stellar mass of 
2.3 x 10 6 Mq, is cored. Common wisdom states that 
one need first wait 10 years for a 30 meter class tele¬ 
scope, such as the Thirty Meter Telescope (TMT), 
to observe such a galaxy. Then one measures the po¬ 
sitions of the stars, waits 5 more years and observes 
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FIG. 1: The precision with which Gaia can measure the 
proper motion of a Sculptor member with Gaia magni¬ 
tude G and color V — Ic = 1-35 


again. This yields proper motions of hundreds or 
thousands stars at a precision of about 4 km/sec, al¬ 
lowing a clean resolution to the cusp/core problem 
around 2030. We will argue that the Gaia satellite 
can yield a 2a hint of the cored/cusped nature of 
Sculptor within 5 years and then, when its position 
measurements are combined with TMT’s first ob¬ 
servations of Sculptor in under 10 years, a definitive 
exclusion of a cusped profile will already be possible. 

In Ref. [5] the authors report the observation of 
1541 objects in the part of the sky occupied by the 
dSph Sculptor, as part of their MMFS survey using 
the Magellan/Clay Telescope [24]. Of these, we con¬ 
sider the 1355 objects for which they assign a 90% 
or greater membership probability in Sculptor. For 
each object they provide the V magnitude and the 
/ magnitude, which we identify with the Johnson- 
Cousins magnitude Ic- We use these to calculate 
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the Gaia magnitude, G , following Ref. [5j 

G = V - 0.0257 - 0.0924(U - I) 

-0.1623(U- I) 2 + 0.0090(R -I) 3 . (1) 

As described in the Gaia Science Performance doc¬ 
ument [7 and updated in [5j to reflect post-launch 
performance, Gaia can measure proper motions with 
an end of mission precision, in p as/yr, of 

= + 680.8(10)°- 4 (G-i5) + 32.7(10)°- 8 (G-i5) 

x0.526 (0.986 + (1 - 0.986)(P - /)). (2) 

There are also position-dependent corrections to this 
precision. Using the location of Sculptor and the 
precision corrections expected [7], we find an im¬ 
provement in this precision of roughly 7 /zas/yr (14 
p as/yr) at G = 17 (G = 18) corresponding to an 
improvement of nearly 20% over the average preci¬ 
sion of 39 fj, as/yr (72 p as/yr). In our calculations 
below we conservatively approximate this improve¬ 
ment to be 15% and so the final precision with which 
proper motions can be measured is 

ct™ = 0.85a™ 0 . (3) 

This function of G is drawn in Fig. [l] for a typical 
color V — I. 

As is explained in Ref. [9], at Sculptor’s ecliptic 
latitude of 36.5° south this improvement is a com¬ 
bination of a geometric parallax factor and a higher 
than average number of transits. However, the num¬ 
ber of focal plane transits is also unusually high for 
Sculptor’s ecliptic latitude, as a result of its longi¬ 
tude and Gaia’s transit pattern. Assuming a dis¬ 
tance to Sculptor of 79 kpc m, the resulting proper 
motion precisions are summarized in the top panel 
of Fig. 0 

We will also determine the precision with which 
the proper motion can be obtained by Gaia including 
a survey of Sculptor by TMT in 2022 out to a radius 
of 700 pc using 10 second exposures of IRIS. This 
will allow an astrometric precision of better than 20 
pas down to magnitude Kab = 20, which includes 
the stars whose positions are well-measured by Gaia. 
IRIS has recently be redesigned to have a 34 x 34 
field of view, so the required observing time will be 
25 hours. In our analysis below, we make the conser¬ 
vative assumption that the astrometric precision at 
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FIG. 2: Number of stars in Sculptor whose proper mo¬ 
tion can be measured with a given precision by Gaia 
alone (top) and Gaia with TMT (bottom). The inclu¬ 
sion of TMT improves the astrometric precison by about 
a factor of 4. 


TMT will be precisely 20 /zas, whereas Gaia deter¬ 
mines positions with a precision of 1.41cr™ p as [7]. 
The combination of Gaia and TMT separated by 6 
years therefore determines the proper motion with a 
precision of 


TMT 


^400 + 2(cr™) 2 ^zas 
6 years 


(4) 


leading to the precisions summarized in the bottom 
panel of Fig. [2] 

One observes that the inclusion of TMT improves 
the astrometric precision by about a factor of 4, an 
approximation which we will adopt in our analysis. 
We have verified that, as the astrometric precisions 
are now smaller than the dispersion, the uncertainty 
with which the velocity dispersion can be measured 
is reasonably insensitive to these precisions. 

In this letter we will answer the following ques¬ 
tions: Given the line of sight velocities of the 1355 
members in the MMFS catalog and proper motions 
measured with the precisions shown in Fig. i (1) To 
what extent can Sculptor’s dark matter density pro¬ 
file be determined? (2) To what extent could the 
CDM paradigm in principle be falsified? 

While there are many approaches to this problem 
in the literature, such as Refs. EHE], our approach 
to these questions will be similar to that of Ref. [4] . 
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However, while we still use an assumed stellar pro¬ 
file to determine the expected dispersion, we per¬ 
form our fit using stars in the MMFS catalog with, 
as described above, proper motions determined with 
a precision given by Eq. ([3]) . The dark matter profile 
p{jr) is assumed to lie in a 5-parameter family, corre¬ 
sponding to the generalization [14] of the Hernquist 
profile [15] 

/ / x b \ ( a ~°)/ b 

' ,(r)= »(r) ( 1+ (^) ) <5) 

where r is the distance from the center of Sculptor. 

We assume both a constant orbital anisotropy 
j3 = 1 — (vg)/(Vr) and a spherical halo. Ref. [2] 
finds that the effect of the first assumption is min¬ 
imal. In Ref. m the authors used the axisymmet- 
ric Jeans equation to determine the elipticity of 6 
nearby dSphs, including Sculptor. They found that 
not only are these systems not spherically symmet¬ 
ric, but in fact their ellipticities are far higher than 
those found in CDM simulations. However, this 
study assumed that the radial and major-axis ve¬ 
locity anisotropies are equal, which in the spherical 
case is equivalent to assuming a vanishing anisotropy 
/3. Indeed, the effect of ellipticity on the line of sight 
velocity dispersion, as described by the authors, is 
quite similar to that obtained by varying /?. Thus it 
seems likely that the statistical significance of these 
high ellipticities would be greatly diminished were 
the restriction on the anisotropies relaxed. This ef¬ 
fect will be compounded by relaxing their other as¬ 
sumptions, such as the parallel orientations of the 
dark matter and stellar distributions. 

Nonetheless, the results of Ref. [IS] demonstrate 
that nonspherical CDM halos do provide a system¬ 
atic modification of the dispersion profiles. This 
shifts the halo parameters obtained using the spher¬ 
ical Jeans equations. Our goal is not to obtain the 
parameters themselves, but merely to determine the 
precision with which they can be determined. By 
assuming that the halos are spherical, we will ob¬ 
tain only a lower bound on the uncertainty on the 
halo parameters. The true uncertainty will be in¬ 
creased by the degeneracy between the spherical halo 
parameters and the anisotropy parameters. If only 
line of sight velocities are available, this degeneracy 
is a serious problem, invalidating such an analysis. 
In a future paper, we will show that with sufficiently 


precise measurements this degeneracy can be broken 
using, in particular, the azimuthal dependence of all 
of the velocity dispersions. However, given Gaia’s 
expected precision we expect that this degeneracy 
will nonetheless increase the uncertainties in the halo 
parameters somewhat beyond those reported below. 

We obtain the radial-dependence of the radial ve¬ 
locity dispersion a 2 = (v 2 (r )) by integrating the 
Jeans equation 

<r 2 (r) = —^3 p s (R)M(R)R^~ 2 dR (6) 

p yr)r j r 


where M(R) is the integrated mass within the radius 
R, approximated to be just the dark matter mass. 
p s [R) is the 3d stellar density at a radius R, not 
just the density of the stars in the catalog, which is 
taken to follow a King profile m with core radius 
7'fc = 0.28 kpc and tidal radius r t = 1-63 kpc [18] . 
The observed line of sight, radial and tangential dis¬ 
persions can be found by integrating ([6]) over the line 
of sight 


°f os (0 

a r( r ) 


<7t( r ) 


P2 D ( r ) 
2 

P2 D ( r ) 
2 

P2D( r ) 



' r 2 \ Rp s {R)a 2 {R)dR 

, P R?) VR 2 - r 2 
' Rp s (RW(R)dR 

, P R 2 ) 


where p\ D (r) is the stellar density integrated along 
the line of sight at r, the distance from the center of 
Sculptor in the transverse plane. 

To determine the precision with which the mass 
profile can be determined, we use the Fisher matrix 


F __ 1 f 1 

2 tr lK(b) + <) 2 99 - de » ) 

( 8 ) 

where the index i runs over the stars in the MMFS 
catalog, the index n runs over the 3 directions los , r 
and t, r, is the projected distance from the center 
of Sculptor to the ith star, 9 a are the 5 parameters 
of the dark matter mass model (|5j) and a™ n is the 
measurement uncertainty on the nth component of 
the velocity of the *th star. The line of sight mea¬ 
surement uncertainties are taken to be 1 km/sec. 
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The derivatives with respect to the parameters 9 a 
are evaluated at the fiducial NFW model [18] 

a = 6=l, c = 3, ro = 0.5, po = 8, /3 = 0 (9) 

where r 0 and p 0 are measured in units of kpc and 
10' Mg/kpc 3 respectively. 

The parameter 0 a can be measured with a pre¬ 
cision By the chain rule, a quantity q 

which depends upon the 9 a may be measured with 
a precision 

W (JM) -&l5- (10) 

We will consider two quantities q. The first is 
M{r i/ 2 ), the mass within the deprojected half-light 
radius rq/ 2 = 375 pc. Using the Jeans equation, this 
can be determined directly from the stellar density 
profile and radial velocity dispersion at r = rq/ 2 , 
their derivatives, and /3 and so is quite robust to dif¬ 
ferent choices of dark matter density profile [HI [20j. 

The second is 

7 (r) = -3 + 47rr 3 ^y (11) 

which is essentially a weighted average of the expo¬ 
nents of the r-dependence of the density at radii be¬ 
tween 0 and r. Note that, in the case of a power law 
p(r), it yields the logarithmic slope. In general it de¬ 
pends on the derivative of M (r), which is determined 
less precisely than M(r) itself, but unlike standard 
definitions it does not directly depend upon the sec¬ 
ond derivative of M(r) and so is nonetheless fairly 
well constrained. Below we will report the fractional 
uncertainty in M(rq/ 2 ) and the total uncertainty in 
7 ^ 1 / 2 )- 

These quantities will be determined in 6 cases. 
The first case just uses the line of sight velocities. 
In the second case we impose the condition c = 3 by 
adding a large number to F cc . This condition is sat¬ 
isfied by both cold and warm dark matter models, 
and in fact by any particulate dark matter model 
with no nongravitational long range interactions. 
Thus, while it is certainly less general, if one is in¬ 
terested in testing CDM then it suffices to impose 
c = 3 and test to see if 7 ~ —1 as is suggested by 
dark matter simulations with the low baryonic con¬ 
tent of the Sculptor dwarf. 

In the next two cases we include proper motion 
measurements from Gaia’s 5 year mission, as has 




67 ( 0 / 2 ) 

LOS only 

69% 

4.1 

LOS only, c = 3 

18% 

1.6 

Gaia 

16% 

0.68 

Gaia, c = 3 

10 % 

0.66 

Gaia+TMT 

10 % 

0.27 

Gaia+TMT, c = 3 

7% 

0.27 


TABLE I: The precision of measurements of the mass 
within ri /2 = 375 pc and the dark matter density slope 
within ri /2 that can be expected with only LOS data 
and also at Gaia with and without TMT. Precisions are 
given using a single component King profile for the stellar 
distribution with and without imposing that c = 3. The 
linearized Fisher matrix approach cannot be trusted for 
entries of order or greater than unity. 

been described above, with and without the condi¬ 
tion c = 3. In the last two cases we add a survey of 
Sculptor by the TMT in 2022. 

These results are all summarized in Table in We 
have also tried including estimates of ro and po from 
CDM simulations, as summarized in Ref. [18]. With 
these constraints we find that Gaia can determine 
7 ( 77 / 2 ) with a precision of 0.57. 

In the standard ACDM paradigm one expects 
mm that a dark-matter dominated galaxy with the 
luminosity of Sculptor will have 7 equal to about -1 
whereas other dark matter models, such as a Bose- 
Einstein condensate or giant monopole model, sug¬ 
gest 7 closer to zero. If 7 indeed is close to zero, then 
by assuming the CDM condition c = 3 and measur¬ 
ing 7^0 one may hope to exclude CDM with 2a of 
confidence when the Gaia mission is completed in 5 
years and with 4 a of confidence when TMT begins 
observing in a decade. One should note however 
that the CDM simulations on which these conclu¬ 
sions rest consider field galaxies, of which Sculptor 
is not an example. 

In Ref. [2TJ it was observed that the Sculptor 
dwarf contains at least two ancient populations of 
stars, a metal rich population near the center and a 
more spatially extended metal poor population. In 
Ref. m the authors found that the stellar density 
and metallicity distribution is well fit by the sum 
of two Plummer profiles of projected half-light radii 
167 pc and 302 pc, having 53% and 47% of the mem¬ 
bers respectively. We have redone our analysis using 
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SM (ri /2 ) /M (j"i/ 2 ) 

< 57 ( 0 / 2 ) 

LOS only 

30% 

1.3 

LOS only, c = 3 

25% 

0.86 

Gaia 

13% 

0.47 

Gaia, c = 3 

13% 

0.44 

Gaia+TMT 

9% 

0.22 

Gaia+TMT, c = 3 

6% 

0.21 


TABLE II: As in Fig. [I] but using a 2-component Plum¬ 
mer model for the stellar density profile. 

this bimodal stellar profile and have found a notable 
improvement in the precision with which 7 can be 
measured, as is summarized in Table [TT] 

If one trusts the equilibrium approximation even 
for fourth moments of the velocities, then these may 
be included without introducing any new free pa¬ 
rameters following the strategy of Ref. [22]. Al¬ 
ternately, as this strategy introduces more observ¬ 
ables without increasing the number of unknowns, 
the proper motion measurements may be used to 
test the consistency of the fourth order equilibrium 
analysis, which may in turn shed light on Sculptor’s 
formation. 

Observation of a core in Sculptor cannot rule out 
WIMPs nor place interesting bounds on their inter¬ 
actions, this would require the observation of a core 
in an ultra faint dwarf (UFD). Gaia cannot precisely 
observe sufficiently faint stars to help TMT with this 
goal. However, in Ref. 23 ] the Hubble Astrometry 
Initiative has been proposed in which Hubble would 
precisely measure the locations of stars in UFDs so 
that future observatories, such as TMT, may deter¬ 
mine their proper motions. This suggests that deep 
observations of the Bootes I UFD by Hubble followed 
by observations by TMT will provide a powerful test 
of the WIMP paradigm. 
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